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Sneutrino-antisneutrino oscillation can be a very useful probe to look for signatures of lepton 
number violation (AL = 2) at the LHC. Here, we discuss the effect of the Lorentz factor 7 and the 
travelling distance L on the probability of the oscillation. We demonstrate that these two parameters 
can significantly alter the probability of the oscillation when the sneutrinos are ultrarelativistic and 
have a very small total decay width. We propose a scenario where these requirements are fulfilled 
and which produces interesting signals at the LHC even for a mass splitting Am as small as 10~ 14 
GeV between the sneutrino mass eigenstates. 
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Oscillation in neutral systems, like K°, B^, B®, 
D°, and u, has been measured and has provided im- 
portant understanding of weak interactions and, in 
the case of neutrinos, revealed their nonzero masses. 
If the neutrino mass is of the Majorana type, it is 
expected that the supersymmetric partners of neu- 
trinos, i.e., sneutrinos, oscillate analogously to the 
neutral meson system. A major difference between 
neutral meson and neutrino oscillations is that, in 
the meson oscillations, similarly to the sneutrino os- 
cillations, the decay of the oscillating particle has to 
be taken into account. In most applications so far, 
the neutral meson can be considered a nonrelativistic 
particle, e.g., this is true in B- factories. The oscil- 
lation of neutrinos is obviously between relativistic 
particles, but the system is qualitatively different 
otherwise, since the neutrinos do not decay. 

Sneutrino-antisneutrino oscillation probes the lep- 
ton number violation (AL = 2) and can be present 
when the neutrinos have nonzero Majorana masses 
[lHlt- It can also provide information on the neu- 
trino sector parameters at the collider environment 
|3, |5|. However, in the derivation of the sneutrino- 
antisneutrino oscillation probability, one usually as- 
sumes that the sneutrinos are produced at rest, as 
in the case of K°-K° oscillation or B°-B° oscilla- 
tion. The situation is different when we produce 
sneutrinos at the LHC energy, and it is not correct 



to assume that they are produced at rest. 

In this paper, we outline the calculation of a 
formula for the sneutrino-antisneutrino oscillation 
probability that is applicable when the sneutrino is 
produced with a very high energy and momentum, 
as is the case, e.g., at the LHC. In the studies of 
sneutrino-antisneutrino oscillation so far, only the 
nonrelativistic case has been considered. We also 
stress the importance of using the correct formula 
in the context of a very interesting supersymmetric 
scenario which can produce spectacular signals at 
the LHC. 

Let us first write down the sneutrino (y) and antis- 
neutrino (z/*) states in terms of the mass eigenstates, 



P> = T^(Pi>+^2», |**> = ;^(|5i>-t|fc». 
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The state \v) at (x,i) becomes 
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The mass eigenstates are \v\) and \v%) with three- 
momenta p\ and P2, respectively. Here, we assume 
that the mass eigenstates move with the same energy 
E but different three- momenta p\ and P2 ■ 

The probability of a \v) oscillating into an \v*) is 



then given by 



p,^ = \(p*mx,t))\ 2 



(3) 



Using Eqs. ([TJ and $Z§, we can expand the proba- 
bility as 



Pc 
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Including the effect of the total decay widths of 
the sneutrino mass eigenstates, one can write down 
the three- momenta pi, with i — 1,2, as pi = 
\fE 2 — m 2 + iTrrii. Here, we assume that the to- 
tal decay widths are the same for the sneutrino and 
the antisneutrino, and the width is denoted by T. 
In addition, m\ and mi are the mass eigenvalues of 
the sneutrino mass eigenstates \v\) and \v-i)-, respec- 
tively. Assuming a very small T and £>T, rri\ , m2, 
we can approximate pi, with i = 1,2, as 

iTrrii 
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(5) 



2{E 2 - m 2 ) 

The last two terms of Eq. (QJ can be expanded us- 
ing standard trigonometric and hyperbolic formulae 
with Im(pi) ~ i^gi and Re(pi) ~ E — ^. Hence, 
we can calculate the probability of a \v) oscillating 
into an \v*) as 
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Here, Am 2 

formula agrees with the formula for neutral meson 

mixing with very large momenta [8(. 

Since the sneutrinos (antisneutrinos) decay, we 
need to look at the integrated probability. Assuming 
mi sa rri2 — to, we get (ip(x, t)\ip(x, t)} ~ e~^~ x and 
f£° dx{ip(x, t)\ip(x,t)) = j^. The integrated prob- 
ability, at a distance L, of a |f/) oscillating into an 
|f5*) is given by 



P(L) 



f Q L dx\(J>*\^(x,t)}\ 2 
Jo°° dx(ip(x,t)\i/j(x,t)) 
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where a = 52S and /3 = ^§-- For a very large L, 
i.e., when La » 1, from Eq. (J7J, we get 

/3 2 
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(8) 
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FIG. 1. Dependence of the integrated sneutrino oscilla- 
tion probability on the Lorentz factor 7 when Al 2 = 1 
and T = 10" 14 GeV. Note that the case of 7 = 1 is 
actually calculated with 7 = 1.0001. 



which is independent of L and where we use the 
relation Am 2 = 2mAm and xo is defined as xo = 
^f 1 [2j. Equation © is the same result as in the 
case when the sneutrinos are produced at rest. Note 
from Eq. (jSJ) that, with La ;^> 1, when xp = 1, the 
oscillation probability P(v — > i>*) is 0.25. On the 
other hand, when x^ ^> 1, P(£ — >• £*) is 0.5. We 
can see from Eq. © that the oscillation probability 
has an exponential suppression factor. 

Next, let us investigate what the effect of the 
Lorentz factor 7 = — on the sneutrino oscillation 
probability is. In order to do this, we must keep the 
length dependence of the oscillation probability for- 
mula (see Eq. (|7|)). Hence, we consider La £0(1). 
Note that the quantity a = ^? = — is the sneutrino 
(antisneutrino) decay width modified by the Lorentz 
factor. In Fig. [TJ we plot the integrated sneutrino 
oscillation probability P(v — > v*) as a function of 
the travelling distance L. The four different lines on 
this plot correspond to four different values of the 
Lorentz factor 7 of the produced sneutrino. We as- 
sume that the total decay width (r) of a 100 GeV 
mass sneutrino is 10 -14 GeV and x„ = 1. It is 
seen from this plot that the oscillation probability 
has a strong dependence on 7 up to a certain value 
of L, and, after that, it saturates and reaches the 
value 0.25, independent of 7 and L. As long as the 
L-dependence is there, for a particular value of L, 
P{y — > v*) is smaller for a larger value of 7. This 
can be understood as follows. Looking at Eq. ©, we 
see that there is a length-dependent exponential sup- 
pression factor that also depends on 7 = E/m. For 



smaller values of 7, this produces a sharper variation 
of the oscillation probability; whereas, for larger val- 
ues of 7, the variation is relatively slow. This is also 
reflected in the variation of the integrated oscillation 
probability with distance, see Fig. [1] 

It is, however, interesting to note that, for a much 
higher value of r(~ lCT 7 GeV) (with x& > 1), the 
value of L is very small («5x 10~ 5 cm) for which 
the oscillation probability saturates (even for 7 = 
50). Hence, for such a large value of T, we can ig- 
nore the effect of 7 or L in the sneutrino oscillation 
probability. 

On the other hand, if the sneutrino (antisneu- 
trino) decay width T is much smaller (i.e., ~ 10~ 14 
GeV or so), the L- and 7-dependences are much 
more pronounced. In such a situation, one should 
use the probability formula given in Eq. ||7J). Such 
small values of the sneutrino decay width are possi- 
ble, for example, in a scenario where the left-handed 
sneutrino NLSP is nearly degenerate to the lighter 
stau LSP and the dominant decay channel for v r is 

*>r^rr+7T + , (9) 

with a total decay width V ~ 10~ 14 GeV. In some 
models with an extra U(1)b-l, the oscillation of 
a right-chiral sneutrino (vr) can be important [9j. 
In such cases, the total decay width of vr can be 
as small as ~ 10~ 14 GeV. The left-chiral sneutrino 
decay width can also be reduced if it has a significant 
mixing with the right-chiral counterpart. 

When the dominant sneutrino decay is v T — > f-j - + 
7T + , one can see a signal pp — > v T Ty — > t^t^ + tt + . 
This produces two heavily ionized charged tracks 
with opposite curvatures when there is no oscilla- 
tion and with same curvatures when there is sneu- 
trino oscillation. We assume that these stau tracks 
can be distinguished from the muon tracks, due to 
the slower velocity of staus. Similarly, one should 
also look at the signal pp — > v*fy — > fy fy + 7r~. 
In this case, when the sneutrino oscillates, one can 
see two same-sign heavily ionized charged tracks due 
to a pair of f-fs. Note that the sneutrino is long- 
lived (decay length approximately a few centime- 
ters), and, hence, one of the staus produced from 
the decay of the sneutrino shows a secondary ver- 
tex which is well separated from the primary vertex. 
This is a very spectacular signal and free from any 
standard model (SM) or supersymmetric (SUSY) 
backgrounds. This parameter region provides dis- 



tinct phenomenology, and one might consider tak- 
ing a point in the region as a benchmark point for 
a general minimal supersymmetric standard model 
(MSSM). 

In order to get an idea about the cross section 
and the branching ratio of the processes discussed, 
we consider a mass spectrum with a v T as the next- 
to- lightest supersymmetric particle (NLSP) and a f\ 
as the lightest supersymmetric particle (LSP). We 
include a tiny i?-parity violating (RPV) coupling 
such that the fi decays outside the detector, leaving 
a heavily ionized charged track. We assume that 
this small RPV coupling does not change the total 
decay width of the sneutrino. 

Nevertheless, below, we list different regions of in- 
terest for the strength of the RPV coupling and the 
width r. In all cases, we assume xo = Am/r xL 1 
(for Xu = 0.2, the maximum oscillation probability 
drops to 2%). 

1. The RPV coupling is very small ( ^S 10~ 8 ) and 
r ^S 10 -13 GeV, as considered in this analysis. 
In this case, the effect of our Eq. (7) is promi- 
nent. 

2. The RPV coupling is larger than what we con- 
sider but does not increase the sneutrino total 
decay width significantly (10 -14 GeV < T < 
10" 13 GeV), and the RPV coupling is - 10~ 7 . 
For such a value of the RPV coupling, the stau 
may decay inside the detector, leaving a heav- 
ily ionizing charged track with a kink. In this 
case, the displaced vertex ( ^ a few mm) from 
the sneutrino will be present. 

3. The sneutrino total decay width is larger 
(but T < 10~ 7 GeV) and the RPV coupling 
is small, ;$ 10~ 8 . In this case, the sneutrino 
oscillation signals remain with charged tracks 
from the long-lived stau. However, in this case, 
the effect of the boost and the displaced vertex 
from the sneutrino will be absent 

4. r < 10~ 7 GeV but the RPV coupling is larger, 

£ 10 -7 . Both the displaced vertex and the 
stau track will be absent, and one has to worry 
about the SM/SUSY backgrounds. 

The mass of the sneutrino is considered to be 



100 GeV and the mass of fi is m-f 



99.7 



GeV. The stau mixing angle is taken to be ir/A. 



Am [GeV] 


10 -14 


io- 13 


10 -io 




Cross section in fb 


Signal 


OS 


ss 


OS 


SS 


OS 


SS 


y/S = 7 TeV 


31.0 


8.1 


20.6 


18.6 


20.3 


18.8 


•/a = 12 TeV 


52.0 


13.6 


34.4 


31.2 


34.1 


31.6 


^/S = 14 TeV 


60.2 


15.8 


39.9 


36.1 


39.4 


36.5 



Am [GeV] 


10 -14 


IO" 13 


10 -io 




Cross section in fb 


Signal 


OS 


SS 


OS 


SS 


OS 


SS 


y/s = 7 TeV 


29.9 


9.3 


19.8 


19.3 


19.6 


19.6 


y/s = 12 TeV 


50.1 


15.6 


33.2 


32.5 


32.8 


32.8 


yfi = 14 TeV 


57.9 


18.0 


38.4 


37.6 


38.0 


38.0 



TABLE I. Cross sections for the OS and SS stau signals 
with several center of mass energies and Am. Here, L = 
0.10 m. The cuts used are mentioned in the text. 



TABLE II. Cross sections for the OS and SS stau signals 
with several center of mass energies and Am. Here, L — 
0.30 m. The cuts used are mentioned in the text. 



The other relevant parameter choices are M\ = 120 
GeV, M 2 = 240 GeV, fj, = -250 GeV, tan/3 = 6, 
m A ° = 600 GeV and A T = 250 GeV. Here, M x 
and M 2 are the U(l) and SU(2) gaugino mass pa- 
rameters, respectively, /i is the superpotential /Li- 
parameter, m A o is the pseudoscalar Higgs boson 
mass and A T is the trilinear scalar coupling of the 
staus. With these values of parameters, the total 
decay width of the sneutrino is T « 1 x 10 -14 GeV, 
while the branching ratio of the decay v r — > f£~ + 7r + 
is 93%. In fact, the branching ratio is greater than 
90% when the mass splitting between the v T and the 
fi is in the range 200-350 MeV. Let us then con- 
sider the production cross sections at the LHC. We 
get the opposite-sign (OS) stau signal pp — > f* f-f 
from both v t t/ and v*f^ productions with an ef- 
fective survival probability (1 — P c «)- The same-sign 
(SS) stau signal pp — > f^ f^ or T-^f-j" we get from 
either v t t± or v*t/ productions with the effective 
oscillation probability (P e ff). 

We select the signal events with the following 
criteria: 1) the pseudorapidities of the staus must 
be \r\ Tl | < 2.5, 2) the isolation variable AR = 
^/(At/) 2 + (Acj)) 2 should satisfy AR > 0.7 for the 
two staus, 3) the transverse momentum of both staus 
must satisfy pT/ > 20 GeV and 4) the /3j should be 
0.3 < f3j < 2.0. The upper limit of f3j reduces 
the muon background considerably. Applying these 
cuts, the cross sections with different center of mass 
energies and different Am are presented in Table U 
for L — 0.10 m. From Table U it is clear that, for 
Am > 10~ 13 GeV, the cross sections almost satu- 
rate. Even putting Am to its maximum value, 10~ 7 
GeV (see Eq. (8) of Ref. [2j), does not change the 
cross sections from Am = 10~ 10 GeV values. On the 
other hand, we can probe down to Am = 10 -14 GeV 
and measure several SS events even with 10 fb -1 lu- 
minosity. Using the SS and OS cross sections, we 



define the asymmetry A = j[ss]+<t(os^ ' ^ or ^ e P ur " 
pose of illustration, we show one value of this asym- 
metry, A = —0.038 ± 0.011, obtained by using ass 
and cros for ^ = U TeV with Am = 10" 10 GeV 
from Table U and assuming an integrated luminosity 
of 100 fb" . This asymmetry A gives direct infor- 
mation about the oscillation probability and is in- 
dependent of initial state parton densities and other 
unccrtainities arising from higher order corrections. 
It is easy to check that P e g = (1 + A)/ 2. By mea- 
suring the value of A, one can calculate the effec- 
tive oscillation probability. For our example, we get 
P cf f = 0.48. 

In Table [TTJ there are the cross sections with dif- 
ferent center of mass energies and different Am for 
L = 0.30 m. These L = 0.30 m values already cor- 
respond to the nonrelativistic oscillation probability 
(i.e., we have Eq. (jHJ at hand). This means that, 
for example, the SS values for y/s = 14 TeV in Ta- 
ble U become ~ 4 — 14% higher if the nonrelativistic 
formula is used. 

If one can measure the three-momentum (\p\) of 
the stau track and the corresponding fj'y at the 
LHC, then one can get an estimate of the stau 
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Ti - 15-, 



10 



11]. The plot of the measured 
stau mass coming from the SS with y/s = 14 TeV, 
Am — 10~ 14 GeV, and L = 0.10 m is shown in Fig. 
[2] All the cuts mentioned in an earlier paragraph 
are used here. The stau momentum and the veloci- 
ties are smeared according to the formulae given in 
Ref. llj . The mass of the decaying sneutrino can 
be measured from the transverse mass distribution 
of the sneutrino. 

In conclusion, sneutrino oscillation is a very im- 
portant tool to look for lepton number violation at 
the LHC. However, at the LHC, the sneutrino can be 
ultrarelativistic, and one should appropriately take 
into account the Lorentz factor 7 = — and the L- 
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FIG. 2. The measured stau mass from the SS with y/s — 
14 TeV, Am = 10~ 14 GeV, and L = 0.10 m. The cuts 
used are mentioned in the text. 



dependence while calculating the probability of os- 
cillation. We have seen that the effect is more pro- 
nounced when the total decay width of the sneu- 
trino is very small (~ 10 -14 GeV), and this can be 
realized in many different SUSY scenarios. A very 
interesting signal at the LHC could be two same- 
sign heavily ionized charged tracks and a soft pion, 
which can probe a mass splitting all the way down 
to ~ 10~ 14 GeV with an integrated luminosity as 
low as 10 ftT 1 for y/s = 14 TeV. In fact, for the 
same mass splitting, it is very evident from Tables I 
and II that, even for y/s = 7 TeV with an integrated 
luminosity as low as 0.5-1 fb _ , one would expect 
to see 4-8 sneutrino oscillation events. 
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